Computational modeling and synthesis of pyridine variants of benzoyl-phenoxy-acetamide with high glioblastoma cytotoxicity and brain tumor penetration

Glioblastomas are highly aggressive brain tumors for which therapeutic options are very limited. In a quest for new anti-glioblastoma drugs, we focused on specific structural modifications to the benzoyl-phenoxy-acetamide (BPA) structure present in a common lipid-lowering drug, fenofibrate, and in our first prototype glioblastoma drug, PP1. Here, we propose extensive computational analyses to improve the selection of the most effective glioblastoma drug candidates. Initially, over 100 structural BPA variations were analyzed and their physicochemical properties, such as water solubility (− logS), calculated partition coefficient (ClogP), probability for BBB crossing (BBB_SCORE), probability for CNS penetration (CNS-MPO) and calculated cardiotoxicity (hERG), were evaluated. This integrated approach allowed us to select pyridine variants of BPA that show improved BBB penetration, water solubility, and low cardiotoxicity. Herein the top 24 compounds were synthesized and analyzed in cell culture. Six of them demonstrated glioblastoma toxicity with IC50 ranging from 0.59 to 3.24 µM. Importantly, one of the compounds, HR68, accumulated in the brain tumor tissue at 3.7 ± 0.5 µM, which exceeds its glioblastoma IC50 (1.17 µM) by over threefold.

Glioblastomas are the most aggressive brain neoplasms with a dismally low 5 year patient survival rate of below 5% 1 . According to the World Health Organization (WHO), glial tumors are classified as grade I and grade II (low-grade gliomas), grade III (anaplastic), and grade IV (glioblastoma) 2 . Current standard-of-care therapies include maximal surgical resection, followed by radiotherapy plus concomitant and maintenance treatment with temozolomide (TMZ) 3 . In addition, a large variety of different genetic and epigenetic modifications have been found in glioblastomas, among which p53, EGFR, PTEN, and IDH mutations are the most common [4][5][6][7][8][9] . However, these validated molecular targets, as well as immunotherapies, including immune checkpoint inhibitors 10 , tumor vaccines 11 , and chimeric antigen receptor T cell (CAR T) therapies 12 , have all been extensively studied, but failed to significantly improve the therapeutic outcome in glioblastoma patients.
There are several reasons why it is difficult to develop more effective glioblastoma therapies: (1) Glioblastomas are characterized by many dysregulated pathways that cannot be blocked simultaneously via a single therapy 13 : (2) Glioblastomas are highly infiltrating and heterogenous tumors that are very difficult to remove by surgical resection without compromising the function of surrounding brain areas 14 ; (3) It is difficult to diagnose glioblastomas in their early stages, therefore, large highly infiltrating and vascularized tumors are often present at diagnosis 15 ; (4) Use of rodent syngeneic and patient-derived models are common to optimize clinical protocols.

Results and discussion
In early stages of drug design it is important to evaluate relevant physicochemical properties of prospective drug candidates. For over two decades, Lipinski's rule of five was the gold standard in drug design 25 . For drugs associated with the CNS, their ability to cross the brain blood barrier (BBB) is the most important characteristic and is not directly addressed by the Lipinski's rule of five. Therefore, two computational scoring systems have been introduced to evaluate the probability of new drug candidates penetrating the CNS. One is CNS Multiparameter Optimization (CNS-MPO) algorithm 26 , which uses 6 physicochemical properties [ClogP (calculated partition coefficient-lipophilicity), ClogD (calculated distribution coefficient at physiological pH (7.4)-lipophilicity), MW (molecular weight), TPSA (topological polar surface area), HBD (hydrogen bond donor at pH = 7), and pKa (-log acid dissociation constant)] for estimating the probability of entering the CNS. The CNS-MPO score range is between 0 and 6, and values ≥ 4.0 have been used as cut-offs for compounds with increased ability to penetrate the CNS 26 . Another scoring system is the blood-brain barrier score (BBB_SCORE) 27 . BBB_SCORE is based on five physicochemical parameters including number of aromatic rings, heavy atoms, MWHBN (a value comprising molecular weight, hydrogen bond donor, and hydrogen bond acceptor), pKa and topological polar surface area. Similar to CNS-MPO, BBB_SCORE also considers the value of 4.0 as the cut-off for acceptable BBB penetration 27 . Consequently, CNS-MPO and BBB_SCORE both serve as complementary algorithms in which CNS-MPO provides information regarding probability of the compound to be found in the CNS, and BBB_SCORE is oriented towards specific physiochemical properties that increase probability of the compound crossing the BBB.
In regard to toxicity, one of the most common hurdles in drug testing is cardiotoxicity triggered by inhibition of the cardiac potassium channel that is coded by the human ether-à-go-go-related gene (hERG). This toxicity test became a mandatory requirement for drug design and development, and can be calculated by the hERG algorithm 28,29 . Based on the compound structure, the hERG algorithm gives a score which is indictive of the compound's inhibitory effect towards hERG. In other words, for the drug to be considered safe, it's hERG IC50 should be significantly higher than its therapeutic IC50 28 .
In order to determine glioblastoma-specific and therapeutically relevant values for all proposed algorithms/ scores, including hERG, we have selected fifteen drugs (see Supplementary Materials, pages 57-79) that are either currently used as glioblastoma drugs, or are in clinical trials with glioblastoma patients 30 . Results in Fig. 1B show that six of these compounds (Fig. 1A) have relatively high water solubility (LogS < 0), low lipophilicity (logD < 3), acceptable values for CNS-MPO and BBB_SCORE (ranging between 3 and 5), and low cardiac toxicity (hERG ≤ 5.5). We used these critical values as a guideline for designing new anti-glioblastoma drug candidates, which are based on the BPA chemical structure present in a common lipid lowering drug, fenofibrate 21,22,31,32 , and in our first glioblastoma prototype drug, PP1 23,24 .
In the search for more effective benzoyl-phenoxy-acetamide (BPA) derivatives, we have decided to explore pyridine variants (Fig. 2). This is because the pyridine moieties are part of a diverse group of compounds with broad pharmacological applications 33,34 . In addition, it is well documented that pyridine-based compounds have a high probability of penetrating the CNS 35 . Herein, we focused on developing and testing five pyridinebase BPA variants, which are grouped according to their structural relationship with pyridine and amide moieties: (I) one with a directly attached unsubstituted pyridine moiety (pyridine-BPA), (II) and one with a     Fig. 3 have low cardiotoxicity (all predicted hERG values are below 5.5), and acceptable water solubility (estimated -logS below 7.5). In addition, CNS-MPO scores for almost all BPA-pyrimidines are over 3, which indicates probability of brain penetration above 50%. Furthermore, these compounds have a BBB_SCORE between 4 and 5, which suggests high probability of crossing the BBB. After determining that the proposed 24 BPA-pyridines have acceptable brain penetration ability, water solubility and low cardiotoxicity, we developed specific preparation methods for each compound. From our previous studies, we demonstrated that the carboxylic group of fenofibric acid (FFA) and its derivatives have exceptionally low reactivity toward nucleophilic acyl substitution. The obvious synthetic pathway for preparation of these compounds is by coupling FFA with corresponding amines 37 . There are two major problems in performing FFA coupling with aminopyridines: (a) FFA has exceptionally low reactivity due to steric constraints generated by the two methyl groups located alpha to the carbonyl group, and (b) amino groups of aminopyridines are exceptionally weak nucleophile due to amino group electron delocalization through the pyridine ring. The difficulty of performing this kind of coupling reaction was recently examined by others 38 . However, FFA can be converted into corresponding fenofibric chloride (Fig. 4), which is sufficiently reactive to aliphatic amines and activated anilines. But, coupling reactions with aminopyridines such as 4-aminopyridine is challenging at best. This is a major drawback to drug design because there are wide applications for pyridine scaffolds in medicinal chemistry 34 . For these reasons, we have decided to explore nucleophilicity of aminopyridines through their computational data.
In molecules of similar structural framework, it is possible to compare frontier molecular orbitals to determine order of nucleophilic and electrophilic reactivity 39 . We have used density functional theory to compute HOMO energies. Results provided in Supplementary Materials (page 153) indicate that aminopyridines should be less reactive (low HOMO) in comparison to aniline while aminophenol is more nucleophilic. Based on our calculations, the coupling of sterically hindered FFA through the reactive fenofibric chloride (FFC) with 4-aminophenol should give a product with a highly isolated yield. Our previous studies 23,40 demonstrated that FFC is a viable, reactive FFA intermediate for amine coupling. For these reasons, we have developed a new synthetic methodology, which is shown in Fig. 4.
There are two distinctive groups of HR-compounds from the viewpoint of amine reactivity and their preparation (a) one with an aromatic heterocyclic ring directly attached to the nitrogen atom (HR66-HR70, HR82, HR83, HR88-HR90), and (b) one with separation of an aromatic heterocyclic ring by one or two aliphatic carbons (HR71-HR81, HR84-HR87). Preparation of the second group of compounds (with reactive amines) is straightforward by mixing FFC with corresponding amines in presence of a base. Because the corresponding amines are more nucleophilic than water and reactions are completed in several minutes at room temperature, these reactions can be performed with an environmentally friendly base such as sodium carbonate in water. Isolated yields are nearly quantitative and this method (Method A) is well suited for small (milligrams) and large  www.nature.com/scientificreports/ (hundreds of grams) scale synthesis. Also, this synthetic procedure does not require any special precautions or preparations. Preparation of the first group, with the heterocyclic aromatic ring directly attached to the amino group, requires certain precautions. Because water is a better nucleophile than these amines' , reaction must be performed in water-free (dry) conditions. In addition, less nucleophilic amines are also prone to oxidation in basic conditions 40 . Because of this, such reactions should be carried out in an oxygen free atmosphere (Method B). Under dry conditions, a pyridine solution of a corresponding aminopyridine was mixed with anhydrous sodium carbonate and kept under nitrogen atmosphere overnight. A separate dichloromethane solution of FFC was also prepared in dry conditions under nitrogen atmosphere. The dichloromethane solution was slowly added to cold (~ 0-5 °C) pyridine suspension of a corresponding aminopyridine and sodium carbonate under a nitrogen atmosphere. Nearly quantitative yields were obtained with reactions performed at 0-5 °C for 3 h, then at room temperature overnight, and finally, at 60 °C for additional 3 h. The product was isolated after solvent evaporation and water addition by simple filtration and extensive water washing. This method gave products in high yield, more the 97% purity, and did not require additional purification by either extraction or chromatography. It is also applicable to milligram and multigram preparation scales for amines with wide ranges of reactivity (Fig. 4). The first group of pyridine-BPA variants are presented in Fig. 5 in which glioblastoma cell viability (CV-based on MTT assay), estimated minimal projection area (MPA), lipophilicity (ClogD), molecular polarizability (PL), as well as, energy of Lowest Unoccupied Molecular Orbital (E LUMO ) and energy of Highest Occupied Molecular Orbital (E HOMO ) were determined. We have included MPA based on studies showing that if a compound does not interact with cell membranes and has a MPA lower than 60 Å 2 , it should be able to penetrate the CNS via passive diffusion [41][42][43] . Therefore, MPA is considered as a better parameter than molecular weight in discriminating the compound ability of entering the CNS 44 .
In the context of comparing compounds with high structural similarity, frontier orbital energies (HOMO and LUMO) can be used in determining which of the similar compounds may have a stronger protein binding ability, and therefore a higher chance of BBB penetration 45,46 .
In this regard, HR67 and HR68 (Fig. 5) both have low HOMO and LUMO energies, along with acceptable lipophilicity (ClogD), and very promising glioblastoma IC50 values (0.59 and 1.17 µM, respectively), indicating that these two compound could become leading candidates for the brain tumor drug development.
The next group of compounds belong to methylenepyridine-BPA (Fig. 6). All of these compounds show promising anti-glioblastoma activity at 25 µM, and their estimated distributions (ClogD) are relatively high, indicating high lipophilicity. MPA values are below 60, suggesting that there is no obstacle to CNS penetration regarding molecular size. Polarizability (PL) is between 40 and 50 indicating that these molecules can adapt to the binding area of a biomolecule through complementary polarization 47 . In the consideration of computed frontier orbital energies (LUMO and HOMO), HR74 should have the best binding ability. In addition, all HR compounds in Fig. 6 are highly cytotoxic at 25 µM except HR75, and corresponding IC50 values for the most promising compounds in this group, HR73 and HR76, are 3.24 and 2.87 µM, respectively.
Next, we asked if the compound activity would change by adding ethylene linker to ethylenepyridine-BPA or by increasing the pyridine's molecular delocalization in benzopyridine-BPA (Fig. 7). As expected, by increasing the number of carbon atoms, either by adding methylene group or additional aromatic ring, lipophilicity increased noticeably (for instance, ClogD for HR83 is 7.29). As the size of the molecule increases, MPA values follow, suggesting that HR81, HR82, and HR83 may have a low probabilities of penetrating the CNS. Considering the energies of both frontier molecular orbitals from three similar compounds (HR78, HR79,  www.nature.com/scientificreports/ and HR80) HR80 should be the most active. Indeed, obtained cell viability (CV) data for these compounds correlated with computational data and suggest that best drug candidate from this group is HR80. Although both benzopyridine-BPAs (HR82 and HR83) have encouraging IC50 values, 1,4 and 2.75 µM, respectively, their computed physical properties such as ClogD and MPA suggest that there is very low probability for these compounds to penetrate the CNS, further supporting HR80 as the best glioblastoma drug candidate in this group. In addition, HR81, which has low glioblastoma-specific cytotoxicity (Fig. 7A) was used here as a negative control for IC50 calculations (Fig. 7C).
We have also explored the importance that chirality may have on the activity of these compounds ( Figure 8). For instance, HR84 and HR85 are structural isomers of HR78 that show moderate cytotoxicity at 25 µM (CV = 26.33). There are noticeable differences between the two stereoisomers: R isomer HR84 being three times more cytotoxic than S isomer HR85 (Figure 8). On the other hand, there is no difference between racemic HR86 (both R&S) and optically pure R isomer HR87. This finding is reasonable because of the pyridine nitrogen proximity to the chiral center. One could argue that pyridine nitrogen binding is sterically diminished in HR85 in comparison with HR84. Because of the position of the pyridine nitrogen regarding chiral center, steric difference is diminished in HR87, which is associated with improved cytotoxicity. The hydoxypyridine-BPAs can exist in both the hydroxypyridine and the amide form with the amide form being preferable by 8   Of the three hydroxypyridine-BPAs, only computed data for HR89 and HR90 can be compared because they are in the amide form while HR88 is in the hydroxypyridine form. All computed parameters suggest that both HR89 and HR90 should easily penetrate the CNS. However, computed frontier orbital energies for HR90 are lower, indicating it should be more effective in penetrating the BBB. In addition, HR90-induced glioblastoma cytotoxicity is nearly 10-fold greater in comparison to HR89, making HR90 the better drug candidate. In conclusion, HR87 and HR90 are the best drug candidates in this group with corresponding, glioblastoma-specific IC50 values of 5.38 and 2.05 µM, respectively (Fig. 8C).
When pursuing new drug candidates, it is very important to envision the possible metabolites and estimate their physicochemical properties, including toxicity 48 . We used computational methods to generate metabolites of all the studied compounds then evaluated their toxicity, solubility, lipophilicity and CNS penetration using the BioTransformer method 49    www.nature.com/scientificreports/ Two pyrimidine variants of BPA, HR67 (PP23) and HR68 (PP21), were subsequently tested for their ability to penetrate artificial BBB model membranes (Fig. 9A,B). In this experiment, BBB permeability values (P) were compared between the two experimental drugs, HR67 and HR68, our prototype drug candidate, PP1 24 , and a positive control (caffeine). We also compared the negative control (FF), for which we have previously demonstrated an inability to accumulate in the brain tumor tissue 22 . Although CNS-MPO scores for HR67 (3.71) and HR68 (3.71) are slightly lower, compared to the CNS-MPO of our prototype drug, PP1 (CNS-MPO = 3.9) 23,24 , these two compounds can cross the BBB model membrane (Fig. 9B). Importantly, we have detected HR68 in the brain tumor tissue at concentrations over 3-fold higher than glioblastoma-specific IC50 for this compound (1.17 µM) (Fig. 9C), further indicating its potential as new glioblastoma drug candidate.
In conclusion, introducing pyridine moieties to the BPA scaffold improves chemo-pharmacological properties of new drug candidates. In particular, water solubility and predicted CNS penetration of these pyridine-derivatives of BPA are higher than in previously studied alkyl and phenolic derivatives of BPA 23,52 . It was also demonstrated here that a properly positioned pyridine moiety in respect to the BPA scaffold increased anti-glioblastoma potency of these compounds with glioblastoma-specific IC50 values being close to 1 µM. Importantly, these specific modifications, which increased molecular flexibility and improved water solubility of the compounds, was achieved without compromising glioblastoma specific cytotoxicity. In addition, stereochemistry of the chiral center close to the pyridine moiety is important because of a difference in three-dimensional orientation of pyridine nitrogen, which can change compound interaction with targeted biomolecules.     , and DMF (few drops) were stirred at room temperature overnight. After approximately 1.5 h, the reaction mixture became light brown. The majority of solvent was removed by distillation at atmospheric pressure and the remaining solvent was removed under Argon flow at room temperature. The resulting solid material was dissolved in dichloromethane (150 ml), under nitrogen atmosphere and with slow stirring, was added to the previously prepared pyridine suspension of 3-aminopyridine and sodium carbonate cooled with ice-water. The resulting suspension was stirred at 0-5 °C for 3 h, then at room temperature overnight, and was followed by stirring at 60 °C for additional 3 h. The resulting solvent was removed under reduced pressure to separate the solid residue. This residue was mixed with water (1 L) and stirred via sonication for 4 h. The insoluble white crystalline product was separated by filtration, extensively washed with water (20 × 100 ml), and dried at 60 °C, under vacuum. The isolated yield was 90% (53.3 g). 1  Trans-endothelial electric resistance (TEER) was measured using a EVOM 2 meter with a STX3 electrode (World Precision Instruments). (B) BBB permeability (P) for the selected compounds was calculated using P = V A ⋅C A /(t⋅S⋅C L ) equation 51 and normalized by TEER coefficient. Data represent average values from 2 independent experiments in triplicates (n = 6) with standard deviation SD. * indicates values significantly different from fenofibrate (negative control), and caffeine was used as positive control. Panel C: HR68 (PP21) tissue concentration evaluated in Foxn1 nude mice bearing intracranial glioblastoma (GBM12). Mice were treated intraperitoneally (ip) with HR68 diluted in 20% cyclodextrin at 15 mg/kg/day and the levels of HR68 in the blood, heart, liver, kidney, spleen, lung, brain and in brain tumor (BT), were evaluated by HPLC, as we previously reported 21,24 . Data represent average values with standard deviation (n = 3). Please note that average glioblastoma IC50 for HR68 is 1.17 µM, and we detected 3.7 ± 0.5 µM of HR68 in the brain tumor tissue.

Method B. General preparation of HR71-HR81 and HR84-HR87 in large scale with sodium carbonate in water as base. Preparation of 2-(4-(4-chlorobenzoyl)phenoxy)-2-methyl-N-(2-(pyridin-4-yl)ethyl)propenamide (HR80).
Fenofibric chloride (FFC) (0.3 mol) was prepared by the following procedure described in Method A from fenofibric acid (95.6 g; 0.3 mol) and oxalyl chloride (63.5 g; 43 ml) in dichloromethane (1 L). Prepared FFC (0.3 mmol) was slowly added dropwise over a period of 45 min at ice-water bath temperature into a magnetically stirred mixture of sodium carbonate (106 g; 1 mol) in water (500 ml) and 1-(pyridine-4-yl) ethanamine (24.4 g; 0.2 mmol) in tetrahydrofuran (500 ml). After addition was complete, the resulting reaction mixture was stirred at room temperature overnight. The reaction mixture volume was reduced by 75% via solvent evaporation under air flow (produced by air pump). The resulting white suspension was mixed with water (500 ml) and the insoluble product was separated by filtration, washed with water (10 × 50 ml), and dried at 50 °C overnight. The isolated yield was 97% (82 g) of pure product. The filtrate was acidified with hydrochloric acid to pH = 2. The resulting white solid precipitate was separated by filtration, washed with water (10 × 20 ml) and dried at 50 °C overnight to give 34.2 g (95% recovery) of fenofibric acid. 1   The resulting mixture was stirred at room temperature in the nitrogen atmosphere for 3 h, followed by stirring under the nitrogen atmosphere at 60 °C for an additional 3 h. After cooling to room temperature, the solvent was evaporated under air flow (produced by air pump) at room temperature yielding a solid residue. This solid residue was mixed with dichloromethane (30 ml) and? water (100 ml). The water layer was discarded, and the organic layer was washed with water (3 × 100 ml), 5% sodium carbonate (3 × 100 ml), and dried over anhydrous sodium carbonate. The drying material was separated by filtration. The volume of the filtrate was reduced to ~ 2 ml, then hexanes (~ 10 ml) were added. The resulting solution was left uncovered at room temperature, allowing the solvent to slowly evaporate. The resulting white crystalline product was separated by filtration, washed with hexane (3 × 3 ml) and air-dried overnight. Isolated yield = 93% (95 mg). Cell culture and viability assays. Human glioblastoma LN-229 cells (ATCC CRL-2611) were maintained as a semi-confluent monolayer culture in DMEM with 1 g/L glucose, sodium pyruvate and L-glutamine (Corning), supplemented with 10% heat-inactivated FBS (Gibco) and P/S (50 units/mL of penicillin and 50 µg/ mL of streptomycin) at 37 °C in a 5% CO 2 atmosphere. Prior to treatment with HR compounds, cells were plated in 96-well plates (BD Falcon) at an initial density of 2 × 10 4 cells/cm 2 . Twenty-four hours after plating, stock solutions of HR compounds were prepared in DMSO, diluted in cell culture medium and added to previously plated cells in triplicate for every experimental condition (final concentration 25 µM). DMSO (0.5%) was used as vehicle control. After 72h incubation, an MTT assay was performed to measure cell metabolic activity (surrogate for cell viability). Following a 1.5 h incubation with 0.5 mg/ml MTT in serum free low glucose DMEM, the resulting formazan crystals were dissolved in 5mM HCl in isopropanol and the absorbance read at 540 nm. Data represent mean values expressed as the percentage of vehicle control ± SD. Phase contrast images of treated cells were taken 72 h after treatment with HR compounds using a BZ-X800 fluorescence microscope (Keyence) equipped with a 20 × objective. The drug dose resulting in 50% inhibition of cell metabolic activity (surrogate for cell viability) was measured using MTT assay, at 72-h time point, and half maximal inhibitory concentration (IC50) was calculated using GraphPad Prism 8.
In vitro model of the blood brain barrier (BBB). The BBB was re-created in vitro using a modified protocol provided by Stone et al. 55 . Briefly, 24-well transwell inserts (Falcon, catalog number 353096) were coated with 10 μg/cm 2 of Collagen Type IV (Sigma) for 24 h at 4 °C. Inserts were washed with sterile water and air-dried for 2 h. Next, the inserts were coated with 2 μg/cm 2 poly-l-lysine (ScienCell) for 1 h at 37 °C, then washed twice with sterile H 2 O and air-dried for 2 h. Primary human astrocytes (1.5 × 10 5 ) and 3 × 10 4 primary human pericytes (both ScienCell) were resuspended in 25 μl of astrocyte medium and pericyte medium (ScienCell), respectively, then combined in a 1:1 ratio for 50 μl total volume. Dried, coated inserts were turned upside down such that the basolateral surface was exposed at the top, and 50 μl of the cell mixture was added to the membrane, covered with the plate lid, and incubated for 2 h at 37 °C to allow cell adhesion. Any medium remaining on top of the membrane was carefully removed before returning inserts to their upright position with the apical surface facing upward, as they were placed in a 24-well plate containing 500 μl per well of astrocyte/pericyte medium  56 were added and incubated for 5 h at 37 °C to allow cell attachment, followed by the addition of an extra 250 μl of TIVE medium. Half the volume of the corresponding media in the lower and upper compartment was replaced with fresh media every third day. Ten days after initial plating, trans-endothelial electric resistance (TEER) was measured using a EVOM 2 meter with a STX3 electrode (World Precision Instruments). The ability of selected HR compounds to pass through the in vitro BBB was tested using inserts with effectively reconstructed BBB as confirmed by TEER values 55,57 .
High performance liquid chromatography (HPLC)-based detection of selected HR compounds. Following TEER measurement, the medium from the apical compartment (insert) of the in vitro BBB model (Fig. 9A) was replaced with 350 µL of fresh TIVE medium containing corresponding compounds [HR67 (PP23), HR68 (PP21), both used at 25 µM. In addition, 25 µM fenofibrate (FF), which does not cross the BBB 22 , was used as a negative control, and 50 µM caffeine was used as a positive control 58 . Plates containing the inserts were returned to the incubator (37 °C, 5% CO 2 ), and after 3 h of incubation, conditioned media from the well and insert (Fig. 9A) were collected. The aliquots (100 μl) of the collected samples were subsequently mixed with 100 μl of 100% acetonitrile, centrifuged (16,000 rpm at 4 °C for 10 min) and supernatants collected for HPLC -based detection of HR67 and HR68. HPLC analyses were performed using an UltiMate 3000 system (Thermo Scientific) equipped with an analytical YMCbasic, 3µm, 150 × 4.6 mm column (octyl silane C8; YMC America, Inc.). Isocratic elution of the compounds was performed using a mobile phase composed of solvent A (50 mM acetic acid in dH 2 O) and solvent B (acetonitrile) mixed at predetermined ratios for each compound (Table 1). All separations were carried out with 5μl sample volume at a flow rate of 1 ml/min, at 25°C. The concentration of each compound was calculated using serial dilutions of the known concentration of the compound separated at the same run with experimental and control samples. After separation, integrated areas under the peak were used to prepare calibration curves and to determine concentration of the compounds.
Intracranial glioblastoma and tissue extraction. Foxn1nu female immunodeficient mice at 6 to 8 weeks of age were used in this study (both male and female are similarly affected by glioblastoma). The mice were inoculated with patient-derived glioblastoma, GBM12-TMZ-resistant, which stably express luciferase reporter 24,59 , were kindly provided by Dr. Sarkaria (Cleveland Clinic, Brain Tumor National Resource) and were cultured and propagated according to the recommended protocols 59 . The cells were injected into the striatum region using 5 μl of PBS containing 1 × 10 5 of the tumor cells guided by the stereotactic approach [1.5 mm posterior to Bregma; 1.5 mm lateral to Sagital suture; 3 mm down from surface] as reported in our previous study 24 . The treatment started when the intracranial tumors were well-established (evaluated by the Optical Image System for small animals (Xenogen IVIS CT). HR68 injection solution was prepared from the 50 mM DMSO stock solution diluted in 20% cyclodextrin (2-Hydroxypropyl-β-cyclodextrin) in sterile PBS and delivered intraperitoneally (ip) at 15 mg/kg. Blood, liver, kidneys, spleen, heart, brain and brain tumor were subsequently collected, solid tissues were washed from blood in PBS, and placed on ice for an immediate sample preparation for HPLC analysis (see above). Tissues were prepared for HPLC by mixing 150 μl of sample tissue (~ 120 mg) that had been mixed with 3 volumes of Methanol : H 2 O mix (4:1), well-blended using TissueRuptor II (Qiagen), and centrifuged at 15,000g for 10 min at 4 °C. Supernatants were collected in 1.5 ml Eppendorf tubes and incubated at 95 °C for 3 min. Following flash cool on ice, samples were centrifuged again at 15,000g for 10 min at 4 °C and supernatants were used for HPLC-based measurement.
Statistical analysis. The data were analyzed with a homoscedastic Student t test. Differences between control and experimental groups were considered significant at P values of ≤ 0.05.

Data availability
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request. License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.